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Abstract. We discuss two-dimensional mapping of the near-infrared emission-line intensity distri- 
butions and kinematics of the narrow-line region (NLR) of NGC 4151, obtained with the Gemini 
Near-Infrared Integral Field Spectrograph, with a projected spatial resolution of w 8 pc. The ionized 
gas intensity distribution follows the projected bi-cone morphology observed in previous optical 
narrow-band images, and its kinematics reveal outflows along the bi-cone. We propose a kinematic 
model in which the gas in the NLR moves at a velocity of « 600 kms^^ up to ^ lOOpc from the 
nucleus. A completely distinct morphology and kinematic structure is observed for the molecular 
gas, which avoids the region of the bi-cone and has velocities close to systemic, and is consistent 
with an origin in the galaxy plane. The molecular gas thus traces the AGN feeding, while the ionized 
gas traces its feedback. 
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INTRODUCTION 

NGC 4151 is the nearest (13.3 Mpc) and apparently brightest Seyfert 1 galaxy, harboring 
one of the best studied active galactic nuclei (hereafter AGN). In the optical, the emitting 
gas of the narrow-line region (hereafter NLR) has the morphology of a projected bi- 
cone, as observed in a previous Hubble Space Telescope (HST) [Oiii]A5007 narrow- 
band image. According to previous studies, our line-of-sight is outside but close to 
the edge of the cones [1, 2, 3], which are oriented along position angle (PA) ^ 60°. 
Optical long-slit spectroscopy reveals outflows along the cones, with the approaching 
side to the SW. In the radio there is a linear structure along PA= 77° that is not aligned 
with the bi-cone [4]. We have observed the inner 3 x Sarcsec^ of NGC 4151 on the 
Gemini North telescope with the Near- Infrared Integral Field Spectrograph ( NIFS) 
[5] operating with the ALTAIR adaptive optics module in order to map the NLR gas 
excitation and kinematics with complete two-dimensional coverage in the near-infrared, 
a spectral region less affected by obscuration than the optical. The spectral range covered 
was 0.94-2.51 jim at two-pixel resolving power ^^5200 and angular resolution of 0.12 
arcsec, corresponding to 8 pc at the galaxy. 

In the top-left panel of Fig. 1, we present a K-band image of the central 60 x 60 arcsec^ 
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FIGURE 1. LEFT - Top: K-band image of the central 60° x 60° of NGC 4151 obtained with the Wilham 
Herschel Telescope. The continuous line shows the orientation of the major axis of the galaxy, while the 
dashed line shows the orientation of the bi-cone. The rectangle shows the region covered by the NIFS 
observations. Bottom: HST [O lll]A5007 narrow-band image of the NLR in the field-of-view of the NIFS 
observations. RIGHT - Intensity maps in the inner 3 x5 arcsec^. Top: [Fell] intensity map; middle: H2 
intensity map; bottom: a coronal line intensity map. 



of NGC 4151, where the large-scale bar can be observed at PA=130°. The position of the 
major axis of the galaxy (22°) is almost perpendicular to the bar. The central rectangle 
shows the field of view covered by the NIFS observations. In the bottom- left panel we 
show the HST [O III] A 5007 narrow-band image of the NLR in the field-of-view of the 
NIFS observations. 

NLR INTENSITY DISTRIBUTIONS AND EXCITATION 

We have maped the intensity distributions in 14 emission lines. Three distinct distribu- 
tions are observed, as illustrated in Fig. 1. Most of the ionized gas is extended to ~ 100 pc 
from the nucleus along the region covered by the known bi-conical outflow, as observed 
in the [Fell] intensity map (top right panel of Fig. 1). The molecular gas (H2) avoids 
the region of the bi-cone, extending from silO to 60pc from the nucleus approximately 
along the large scale bar, almost perpendicular to the bi-cone axis (middle right panel 
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FIGURE 2. Position-velocity diagrams of [S III] 0.9533 jUm emission for 0.2 arcsec wide pseudo- 
slits oriented at PA=60° along the bi-cone axis and offset perpendicular to the pseudo-slit by —0.8 
arcsec to -1-0.8 arcsec as indicated in each frame. An average continuum image extracted between 
±(2000 — 2500) kms^' has been subtracted. Crosses mark the location of the nucleus in the central 
pseudo-slit and the "equivalent" location for the other slits. 



of Fig. 1). The molecular gas appears to be destroyed by the intense nuclear radiation 
escaping along the bi-cone, and its intensity distribution suggests an origin in the galaxy 
plane. The coronal line region is only barely resolved (bottom right panel of Fig. 1), 
consistent with an origin in the inner NLR. 

The line ratios [Fe ll]1.257jUm/[P ll]1.189/im and [Fe ll]1.257jUm/Pa/3 oftheNLRof 
NGC4151 correlate with the radio intensity distribution, mapping the effects of shocks 
produced by the radio jet on the NLR. These shocks probably release the Fe locked 
in grains and produce the observed enhancement of the [Fe ll] emission at arcsec 
from the nucleus (top right panel of Fig. 1). In these regions, we obtain electron densities 
He ~ 4000 cm^^ and temperatures Te ~ 15000 K for the [Fe ll]-emitting gas. For the H2- 
emitting gas, we obtain much lower temperatures of Texc ~ 2100K and conclude that 
the gas is in thermal equilibrium. The heating necessary to excite the molecule may be 
due to X-rays escaping through the nuclear torus [8]. The escape of radiation through 
the torus is also required by the observation that the bi-cone does not have a sharp apex, 
but rather show ionized gas emission in all directions around the nucleus. 

GAS KINEMATICS 

In Fig. 2, we present a position-velocity (PV) diagram for the [S ill] emission line along 
the bi-cone axis. We modelled the NLR kinematics as a truncated hollow conical outflow 
at a constant velocity of 600 km s ^ We did not see acceleration or deceleration in 
our data, as proposed by Das et al. [3] from the modeling of the [O ill] emitting- 
gas kinematics. The model PV diagram is shown in Fig. 3. Besides the outflowing 
component, the ionized gas shows two components at velocities close to systemic. One is 
gas from the galaxy plane, ionized by radiation escaping through the torus, as discussed 
above, and the origin of the sequences at zero velocity in Fig. 3. The other is emission at 
zero velocity generated by the interaction between the radio jet and gas from the galaxy 
plane (the three knots in Fig. 3). 

We estimated a mass outflow rate of ^ IMQyr^^ along each cone, which exceeds the 
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FIGURE 3. Simulated position-velocity diagrams for a truncated hollow cone model with constant 
velocity of 600 kms^' outflow. Besides the bi-conical outflow zero velocity components are shown for 
the emission associated with the radio jet (the three knots) and for the galactic disk component (sequence 
at zero velocity). 



inferred black hole accretion rate by a factor of ~ 100. This can be understood if the 
NLR is formed mostly by entrained gas from the circumnuclear interstellar medium by 
an outflow that probably originates in the accretion disk. 

The H2 emission, besides presenting a distinct intensity distribution, also has a distinct 
kinematic structure from the ionized gas. The velocities are close to systemic and 
consistent with an origin in the galaxy plane. This hot molecular gas may only be the 
tracer of a larger reservoir of colder gas which represents the AGN feeding, while the 
ionized gas outflow represents the AGN feedback. 

More details on the NLR excitation and kinematics can be found in Storchi-Bergmann 
et al. [6] and Storchi-Bergmann et al. [7], respectively. The analysis of the nuclear 
continuum - which reveals the presence of an unresolved infrared source that can be 
attributed to a dusty torus — is discussed in Riffel et al. [8]. 
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